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We investigate the effects of wave localization on the delay time r (frequency sensitivity of the 
scattering phase shift) of a wave transmitted through a disordered wave guide. Localization results 
in a separation t — x + x' of the delay time into two independent but equivalent contributions, 
associated to the left and right end of the wave guide. For = 1 propagating modes, x 
x' are identical to half the reflection delay time of each end of the wave guide. In this case the 
distribution function P(r) in an ensemble of random disorder can be obtained analytically. For 
JV > 1 propagating modes the distribution function can be approximated by a simple heuristic 
modification of the single-channel problem. We find a strong correlation between channels with long 
reflection delay times and the dominant transmission channel. 



I. INTRODUCTION 

In this paper we characterize localization of randomly 
scattered waves by means of a dynamical quantity, the 
delay time r. 

Wave localization is perhaps the most striking effect 
of multiple random scattering — in a wave guide ge- 
ometry, it results in the exponential attenuation of the 
transmitted intensity I{L) cx exp(— 2L/^) for lengths L 
of the wave guide greater than the localization length 
f , even in the absence of absorption. Localization was 
first investigated in mesoscopic systems Recently 
the undertaking of its realization and observation for mi- 
crowaves MM and optical waves |ic|] has attracted a lot of 



interest. It is still under debate |11 I4I whether some of 



these observations are due to localization or absorption. 

The delay time r — d4>/du! is the frequency sensitiv- 
ity of a scattering phase shift 0, and has been identified 
by Wigner as a measure of the exploration time of 



the scattering region (see also Refs. [p_4 



L5|). Recent ex- 



periments have succeeded in the direct measurement of 
the so-called single-mode delay time for specified incident 
and detected modes, both for microwaves |jl^ and optical 
waves [|7j. (The attribute 'single- mode' means here that 
only one of the N propagating modes is excited, and only 
one mode is selected for detection, but does not imply 
any restriction of N itself.) These experimental efforts 
have promoted the single-mode delay times to quantities 
of interest in their own right. The measurements have 
been performed with wave guides shorter than the lo- 
calization length, and their outcome can be successfully 
described by diffusion theory Q. That does not mean 
that wave localization is of no interest in this context — 
note that the experiments on localization and delay times 
have been performed on the same sorts of sample, by the 
same groups. 

Theoretical work on the localized regime has mostly 



concentrated on the delay times of the reflected signal 
Jl9|-p6|. Some aspects for the transmission delay time 
problem for a single propagating channel (iV = 1) have 
been studied in Ref. ||27|, where it was found that the 
distribution of r has a universal quadratic tail, P{t) cx 
T~^, for large t. This tail eventually crosses over into a 
log-normal tail, at some large value Tc that increases with 
the system length — even though the tail is irrelevant for 
the direct experimental or numerical investigation of the 
distribution itself, it is reflected in physical properties of 
mesoscopic systems (for a review see Ref. Ref. |2^] 

also addressed the properties of a delay time weighted 
by the transmission coefficient, which is relevant for the 
conductance of mesoscopic wires. 

In this work we investigate the distribution of the 
transmission delay time r in the localized regime. It will 
turn out that the transmission and reflection problem are 
closely related for = 1. The transmission delay time is 
then the mean of the reflection delay times of both ends 
of the wave guide, and the exact form of the limiting 
distribution function P(t) for L — > 00 can be found an- 
alytically. At finite length the result is applicable in the 
range < t < Tc- Because Tc is very large in the localized 
regime, this covers the range of delay times which is rele- 
vant for direct experimental observation and comparison 
with numerical simulations. 

For > 1 there is still only one relevant transmission 
channel. Consequently, once again localization results in 
a separation of the transmission delay time into two inde- 
pendent but equivalent contributions from both ends of 
the wave guide. Moreover, one of the contributions only 
depends on the excitation mode, while the other only de- 
pends on the detection mode. However, the transmission 
delay times are no longer directly related to the reflec- 
tion delay times. Nevertheless it is possible to obtain the 
distribution function of single-mode delay times approx- 
imately by a heuristic modification of the single-channel 
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FIG. 1. Quasi-one dimensional wave guide filled by a dis- 
ordered medium and illuminated by a monochromatic plane 
wave. The scattered wave acquires a scattering phase shift 
(j). We investigate the frequency sensitivity (delay time) 
T = d(j)/diu for the transmitted wave. 



problem. 

Although there is no direct relation to the reflection 
problem for the individual single-mode delay times and 
A'^ > 1, there exists an intensity- weighted combination 
of all delay times which is more closely related to the re- 
flection problem. This combination involves the orthogo- 
nal transformation matrix from the basis of transmission 
channels to the eigenvectors of the Wigner-Smith time- 
delay matrix. From our numerical simulations we find 
a strong correlation of the dominant transmission chan- 
nel and the channel with the largest Wigner-Smith delay 
time. 

The paper is organized as follows. In Section | we 
provide the necessary background material which will be 
used later on in the investigation of the transmission de- 
lay times. This includes a short review of the diffusive 
regime and the reflection delay times in the presence of 
localization. In Section |l| we discuss the case iV = 1 of 
a single-channel wave guide and calculate the distribu- 
tion function of the transmission delay time analytically. 
Section IV is devoted to wave guides with more than 
one propagating channel. We will first discuss the single- 
mode delay times and compare the distribution from a 
numerical simulation with the analytic expression that 
arises from the heuristic approximation. Then we turn 
to the weighted combination of all delay times and use it 
to investigate the relation of the dominant transmission 
channel with the channel associated to the largest delay 
time. 



II. BASIC CONCEPTS 



A. Wave-guide geometry 



Fig. depicts a quasi-one dimensional wave guide 
(length L much larger than the width) which is filled by 
a medium with randomly placed scatterers (mean free 
path I). We assume that there is no absorption and no 
inelastic scattering inside the wave guide, and consider 
a monochromatic scalar wave (disregarding polarization) 
for simplicity. Also we assume that time-reversal sym- 
metry is preserved, as is appropriate for the propagation 
of light in absence of magneto-optical effects. 



The number N of propagating modes at frequency uj 
equals the number of transversal excitations inside the 
wave guide, and is given by A^ = ■nA/X'^ for a wave 
guide with openings of area A (here A = w/c is the wave 
length and c is the propagation velocity of light). In 
the numerical simulations we will work with a planar 
wave guide of width W <^ where A^ = 2W/\. For 
a unified description we introduce the scattering time 
7 = al/c, with the coefficient a — 2 (7r^/4, 8/3) for one- 
dimensional (two-dimensional, three-dimensional) scat- 
tering inside the quasi-one dimensional wave guide, and 
the relative length s = a'L/l, with a' = 1/2 (2/ it, 3/4). 
The localization length is then given by ^ = {N + 1)1 /a'. 



B. Scattering formalism 

The number N of propagating modes inside the wave 
guide corresponds to the number of independent incident 
modes close to each opening of the wave guide. In experi- 
mental practice these modes can be chosen as plane waves 
with discretized propagation direction, and mode selec- 
tion is realized by the choice of the positions of source 
and detector. In such a single-mode experiment, the wave 
guide is probed by external illumination with amplitude 
^I'm in mode m, and the transmitted or reflected signal 
$„m is detected in mode n, with n,m = 1, . . . , 2N . (The 
modes with index n,m = 1, . . . , N are associated with 
the left end of the wave guide, while the remaining modes 
pertain to the right end of the wave guide.) The numbers 

Snm = «'nm/*m (1) 

form the elements of the 2N x 2N scattering matrix 



S = 



r t' 
t r' 



(2) 



with four N X N dimensional blocks which correspond 
to reflection or transmission with the incident radiation 
from the left (r, t) or from the right (r', t'). The scatter- 
ing matrix is unitary due to flux conservation in the ab- 
sence of absorption, and only depends on one frequency 
because there are no inelastic processes. Furthermore, 
the scattering matrix is symmetric due to time-reversal 



symmetry, hence t' = t 



, and r' 



r.iT 



A useful representation of the scattering matrix is the 
polar decomposition 



u 








T 



u 
V 



(3) 



with unitary matrices u and v and the diagonal ma- 
trix T = diag (Ti, . . . , T/y) of transmission eigenvalues 
(eigenvalues of tH). For convenience we order them by 
magnitude, Ti > T2 > . . . > Tn- 
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C. Intensity and delay time 



The elements of the scattering matrix can be written 



as 



^exp(i(? 



(4) 



where Inm is the detected intensity for unit incident in- 
tensity and is the scattering phase shift. The single- 
mode delay time is defined as the derivative of the scat- 
tering phase shift with respect to frequency, 



also be understood from the uniform distribution of the 
matrices u and v in the group of unitary matrices U(A'^). 
For large N, the elements of u and v can be considered 
as random Gaussian numbers with variance {\uim\'^) = 
(|wimp) — l/N, and the Rayleigh distribution (p) follows 
from the central-limit theorem. 

The distribution function of the delay time is given by 
@@ 

P{rnm) = #t[Q + (r„™/(r) - l)^]-^/^. (8) 



(5) 



In transmission (5 = 2/5 and (r) ~ 7*^/3, while in reflec- 



Its interpretation as an exploration time of the medium 
stems from the short-wave-length limit. The phase can 
then be approximated by the classical action 5ci of trajec- 
tories (there may be several) which satisfy the boundary 
conditions of the incident and detected modes. Accord- 
ing to classical mechanics, the derivative dSd/du; of the 
phase with respect to frequency (energy) equals the clas- 
sical propagation time through the medium. 



D. Ballistic case 

In the ballistic regime s <C 1 the wave is transmitted 
without any attenuation, and the modes can be chosen 
easily such that each incident mode m is strictly associ- 
ated with a transmitted mode n'{m), namely, by using 
the reflection symmetry of the wave guide (exchanging 
left and right). The intensity is then given by Inm = Snn' , 
and the delay time is r„m — (5„„/L/cm, where is the 
longitudinal propagation velocity in mode m. The aver- 
age over all modes is (L/cm) — js. 



E. Diffusion theory 

Diffusion theory applies when the length L of the wave 
guide exceeds the mean free path I but is less than the 
localization length ^. The fluctuations of the intensity 
Inm for given m and varying n result in a speckle pat- 
tern of bright and dark spots, which is described by the 
Rayleigh distribution 



P{Inm) = -^cxp (-/„,„/(/)). 



(6) 



The mean intensity per mode is (/) = {T)/N in trans- 
mission and (/) = (1 — {T))/N in reflection, where M 



(r) = l(trttt) = (l 



(7) 



is the mean transmission probability. For the special case 
n = m in reflection the mean intensity doubles due to 
coherent backscattering The speckle pattern can 



tion Q — 2s/5 and (r) 
in reflection, see Ref. |23]). 



275/3 (for ballistic corrections 



F. Localized regime 

In the localized regime i ^ ^ the transmission 
eigenvalues r„ become exponentially small, with well- 
separated, self-averaging exponents — (lnT„)/L — 2n/£_. 
Transmission is dominated by the transmission channel 
with eigenvalue Ti , which is exponentially larger than all 
the other transmission eigenvalues. In terms of the polar 
decomposition (||), 



tr. 



'TiVinUir 



Inm — Ti\vinUln 



(9) 



For large N the complex numbers wi„ and ui„i again can 
be considered as Gaussian random numbers. For fixed 
incident mode m and within a given disorder realization 
(fixed Ti), this results again in the Rayleigh distribution 
(I) for Inm, with (/) = Til vim\ ■ If one also averages 
over the incident mode, however, one finds 



P{,Inni) 



2N' 
~I\ 



Ko{2N^Inm/Ti), 



(10) 



with Kq a modified Bessel function of the second kind. 
This deviates from the Rayleigh law, obviously because 
the central limit theorem no longer holds due to the large 
relative differences between the transmission eigenval- 
ues. The reflected intensities /,„„, however, still follow 
the Rayleigh distribution with (/) — l/N, since they 
are governed by the non-fluctuating reflection eigenval- 
ues i? , = 1 - Ti « 1. 

Because transmission becomes negligible, the reflection 
matrices r = u^u and r' = —v^v become unitary. The 
single-mode delay times of reflection can then be related 
to the Wigner-Smith delay times fi, f/, which are the 
eigenvalues of the Wigner-Smith matrices 

. *dr t /^n T * ^""^ ^ 
q = —ir' -r- = u ' I 2 Im u — — | u, 



duj 
duj 



du) 



i 2Inii;* I V, 

duo 



(11) 
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respectively (for details of the relation refer to Refs. 

The two sets of Wigner-Smith delay times are inde- 
pendent and equivalent. In terms of the rates /ii = 
the joint distribution function is given by the Laguerre 
ensemble E3l 



-7(^+1)^. 



(12) 



where the step function 8(2:) = for x < and Q{x) — 1 
for a; > 1. Eq. ( p^ generalizes earlier results for = 1 
||9|-|2| to arbitrary N. 

We order the delay times by their magnitude, fi > 
f2 > . . . > ttv • Of special interest is the largest delay 
time fi, which is known to dominate the statistics of 
the reflection delay times p^ , p5[ , although to a lesser 
extent than Ti determines the transmitted intensity. Its 
distribution follows from a result by Edelman po|] for the 
smallest /i in the Laguerre ensemble and is given by 



Pin) 



jN{N +1) 



exp[-77V(iV + l)/fi] 



(13) 



The mean (fi) diverges because of the quadratic tail for 
large f 1 . These large fluctuations are a signature of local- 
ization 1^,^,^ 32 1 , and have been interpreted as explo- 
ration of the localized regions deep inside the wave guide. 
Our result for the transmission delay tim e will support 
this interpretation: We will see in Section IV D that the 
corresponding eigenvector of the Wigner-Smith matrix is 
correlated with the dominant transmission channel. 



III. SINGLE-CHANNEL WAVE GUIDE 

The distribution of the transmission delay time T12 for 
a single propagating mode {N = 1) has been investigated 
previously in Ref. where it was found that P(ti2) oc 
r{^^ for large T12. In this Section we will calculate the 
distribution function analytically, for all T12. 

For = 1 the scattering matrix is a 2 x 2 matrix, hence 
the transmission and reflection elements t = uvVt, 
r = — T, and r' — — T reduce to com- 

plex numbers, while the matrices u, v, of the polar de- 
composition are now unimodular complex numbers. The 
single-channel case is special because the transmission 
delay time 

T12 = Imu h Imu — = (14) 

auj duj 2 

is directly related to the reflection delay times 

Til = Imr"^ — = 2ImM"^ — , (15 
dbj dw 

T22 = Imr'-i^ =2Imi;-i4^. (16) 
du du) 



The relation holds for all lengths (it does not require 
localization) , and also can be derived from the condition 
of unitarity of the scattering matrix. 



rt* 



tr' 



: 



duj 



{rt* +tr'*) = 0. 



(17) 



It is convenient, also in view of the case TV > 1 to be 



discussed in Section IV, to introduce the quantities 



X = ImM — — , X = fiii''^ • 
duj duj 



(18) 



In the localized regime, the reflection delay times are de- 
termined by scattering in non-overlapping regions close 
to each end of the wave guide. Hence x ^^'^ x' be- 
come independent, and their joint distribution function 
P{XtX') = P{x)P{x') factorizes. The reflection delay 
times Til = 2x, T22 = 2%' equal the Wigner-Smith delay 
times fi, f{, respectively. The function 
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exp(-7/x)e(x) 



(19) 



[and equivalently P{x')\ hence follows from the Laguerre 
ensemble, Eq. (p^, for iV = 1. The derivation in the 
framework of one-dimensional scaling theory is briefly 
recapitulated in Appendix ^ P{x) eventually is deter- 
mined by the requirement that it becomes independent 
of length in the localized regime, which results in the 
stationarity condition 



7C' 



dP 
'dL 



d_ 

dx 



d . 
dx 



P = 0. 



(20) 



In Section IV B we will propose a slightly modified ver- 
sion of this equation for the case A^ > 1. 

From Eq. ( p^ ) , the distribution of the transmission de- 
lay time T12 = X + x' is then found by integration. 



Piri2) 



dx ■ 



7 



■ exp 



[ti2 - a;)2 
T12 



Ko 



exp[-7/ X- j/ (ti2 - x)] 



27 
T12 



27 
T12 



(21) 



In Fig. ^ this prediction is compared with the result 
of a numerical simulation of random scattering in a pla- 
nar single-channel wave guide. In these simulations the 
Helmholtz equation is solved on a square lattice. In 
terms of the lattice constant a, the width of the wave 
guide is W = 3 a, and the wave length is A = 4 a, giving 
rise to a single propagating mode. Disorder is modeled 
by a random on-site potential, with localization length 
^ = 4Z = 54 a. The scattering rate 7 is determined from 
the ballistic regime. We find perfect agreement between 
Eq. (^) and the numerical simulations, without any free 
parameter. 
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FIG. 2. Distribution of transmission delay time r for a 
single-channel wave guide. The analytic result (^l|) (curve) is 
compared with the results of a numerical simulation of ran- 
dom scattering in a single-channel wave guide. 



IV. MULTI-CHANNEL WAVE GUIDE 

Now we turn to the case > 1 of more than one 
propagating mode in the wave guide. We first show that 
the delay times separate into two independent contribu- 
tions and discuss some consequences. Then we turn to 
the distribution function P{Tnm) and propose an approx- 
imation, based on a heuristic modification of the case 
iV = I, which agrees well with the result of numerical 
simulations. Finally, we investigate the correlation be- 
tween the transmission channel with eigenvalue Ti and 
the eigenvector of the Wigner-Smith matrix with eigen- 
value fi. 



A. Separation rule 

For the transmitted intensity it is sufficient to con- 
sider the reduced form tnm — ^^inWimVTi, Eq. (^, of 
the transmission-matrix elements in the localized regime. 
Under the additional assumption (which is validated by 
the numerical simulations) that the coefhcients Vkn, Ukm, 
k ^ 1, do not depend much more sensitively (by large fac- 
tors oc e^'"/^) on frequency than the elements wi„, Wim, 
this form also can be used for the delay times, which then 
separate into two contributions. 



2 Xm 

= Im- 



1 dui„. 
duj 



x'n = Im 



1 dvin 



(22) 
(23) 



The contribution Xm only depends on the mode index 
m of the incident mode, while x'n only depends on the de- 
tected mode n. This gives rise to strong correlations be- 
tween the delay times for each disorder realization: They 
obey the relations 



+ Tki — ni + . 



(24) 



The dependence on the mode indices suggests that x 
and x' are independent and equivalent, and that they are 



determined by scattering within a couple of localization 
lengths close to the associated opening. This is also sug- 
gested by the fact that Xm only depends on the matrix 
M, while x'n only depends on the matrix v. These matri- 
ces, on the other hand, determine the reflection matrices 
r = u^\/l — T u w u^u and r' w —v'^v, which can be 
considered as independent in the localized regime. (The 
approximation T = corresponds to neglecting the in- 
fluence of the opposite end of the wave guide, which is 
far away). However, that might be deceptive — note that 
although u and v give r and r', they are themselves not 
uniquely determined by r and r' in this approximation: 
E. g., the same reflection matrix r can be obtained from 
ou, with o an arbitrary orthogonal matrix. The matrix u 
only can be determined uniquely from r if we also use 
the information in T, which depends on the opposite 
end of the wave guide. We will demonstrate now that 
X and x' nevertheless become ind epend ent in the local- 
ized regime. However, in Section [V C we will see how 



degrees of freedom similar in nature as o reflect in the 
statistical distribution of the delay times. 

In order to demonstrate that x and x' are indeed inde- 
pendent, we cut the wave guide into two parts (associated 
with subscripts i = 1,2), still requiring that the lengths 
Li ^ 1^. The well-known composition rule 



t = t2{l - r[r2)-Hi 
and the relations 1 ^ Ti^i 3> Tfc^^i.i yield 



(25) 



(26) 
(27) 



Note that Ti is indeed real. This gives r„m = Xm.i+x'n 2; 
i.e., Xm = Xm,i independent on part 2 and x'n = x'n,2 
independent on part 1. 



B. Distribution of delay times 



The considerations in the previous Section IV A also 



show that the statistical distribution of r becomes inde- 
pendent on length ("stationary") for L ^ ^, because the 
distribution P{xm) for length L is identical to P{xm,i) 
for length Li < L, and analogously for x'n- 

The stationary distribution P(t) is plotted in Fig. ||, 
for N = 2 and TV = 30 propagating modes in the numer- 
ical simulations (corresponding to different widths W of 
the wave guide). The distributions collapse onto a sin- 
gle curve when the delay times are rescaled by a factor 
jN{N +1). The distribution is however qualitatively 
different from the result for = 1. Most notably, a 
tail cx also extends into the region of negative de- 
lay times, while the delay times for = 1 are strictly 
positive. 

An analytic treatment of the transmission delay time 
problem for many channels is notoriously difhcult. In the 



5 




FIG. 3. Distribution of transmission delay time r„m 
for multi-channel wave guides with N = 2 (open dots) and 
N = 30 (full dots). The analytic prediction from Eq. (H) 
(curve) with Y = Z = ^N{N + l)/2 is compared with the 
results of a numerical simulation of random scattering in a 
planar wave guide. 

framework of one-dimensional scaling theory, the evolu- 
tion of X couples to all elements of u and du/duj, which 
makes a complete analytic solution impossible. Inspec- 
tion of the complicated full set of evolution equations 
that appear in this approach, however, suggests the fol- 
lowing approximation for the stationarity requirement of 

Pix)-- 



z") Pix) = 0, 



with the solution 



P{x) 



Y exp[(r/Z) arctan(x/^)] 
2(^2 + Z2) sinh(7ry/2Z) ' 



(28) 



(29) 



For Y — J and Z = 0, the stationarity condition re- 
duces to Eq. (|20| ) for = 1. For > 1, the appearance 
of Z can be traced back to the additional degrees of free- 
dom in u^du/duj, especially also to the real part of this 
matrix (the real part vanishes for N — 1). This will be 



furt her di scussed in the following two subsections IV C 
and |IVD| . The factor Y/Z in the exponent of Eq. (|2|) 
determines the asymmetry of the distribution for positive 
and negative values of x- 

The full set of evolution equations suggests that Y ~ 
Z ~ 'yN{N + l)/2, up to numerical factors which cannot 
be derived without solving the original problem. This 
is also the order of magnitude of r„m at the border of 
diffusion and localization, see Eq. (||). In Fig. | we have 
plotted the distribution of r = x + x' which follows from 
Eq. (|2|) for Y = -fN{N + l)/2 and Z = -fN{N + l)/2. 
The comparison with the numerical data shows that the 
numerical factors are close to unity. 



times. Now we discuss to which extent these two prob- 
lems are linked for > 1. 

Due to symmetry, the scattering matrix can always be 
written as 5 = U^U. In terms of the matrices of the 
polar decomposition, we can choose 



U 



(r/2p)i/2w 

-^(p/2)l/2^ 



(W2) 
t{T/2p) 



1/2, 



1/2, 



(30) 



with p — 1 — vl- 
approximated by 

U = 



T . In the localized regime, U can be 



2 ' 

iv 



(31) 



The first index of the matrix U is decorated by the trans- 
mission amplitudes, hence U relates the scattering states 
to the transmission channels (each transmission channel 
is characterized by two vectors: a row of u which con- 
nects it to the scattering states on the left and a row of 
V which connects it to the scattering states on the right). 

The Wigncr-Smith time-delay matrix of the total scat- 
tering matrix is 



auj 



.dU 



duj doj 



(32) 
(33) 



From the unitarity of U is follows that Q' is real and sym- 
metric, and hence diagonalized by an orthogonal matrix, 
which we write in block form 



O = 



On 012 

021 022 



(34) 



In this block form we denote the set of eigenvalues as 
diag(f,f'). The matrix O diagonalizes the Wigner- 
Smith matrix Q in the basis of transmission channels, 
given by [/, and hence relates the transmission channels 
to the eigenvectors of the Wigner-Smith matrix. 

It is consistent to assume that O is almost block- 
diagonal, with off-diagonal elements 012, 021 of order 
\/Ti. From Q' = O diag (f , f ')0-^ we indeed obtain un- 
der this assumption the relations 



oiifoii = 2Imu*— — , 
dui 



.du" 



022t'o22 — 2Imw^ 



dv^ 
dio 



(35) 
(36) 



02lfoJ^ + 022f'0i2 = Re 



dio dio 



(37) 



C. Relation to the reflection problem 

For iV = 1 we could relate the problem of transmission 
delay times directly to the problem of reflection delay 



Comparison with Eq. ( |11D shows that on and 022 diago- 
nalize the Wigner-Smith matrices of the reflection prob- 
lem (Section II F), however, in the special basis of trans- 
mission channels which is not fixed by reflection alone. 
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The matrices 012 and O21 are related to frequency deriva- 
tives of u and v which do not feature in the reflection 
problem at all. Moreover, because they appear as off- 
diagonal elements of O, these matrices connect the co- 
efficients of the transmission channels from one side of 
the wave guide to Wigner-Smith eigenvectors of reflec- 
tion from the other side. 

According to Eqs. ( ^5|) and (|6|), the eigenvalues of Q 
can be approximated by the two sets f, f of Wigner- 
Smith delay times of the reflection matrices. The trans- 
mission block of 



duj 



(38) 



corresponds to 



Re- 



[v (vToii -I- 2io2i)fo{iu\ 



+ Re 



[v'^ 022f' {O22VT + 2iof2)u]r 



(39) 



Note that the separation ( p3| ) of the transmission delay 
time into two contributions which only depend on the 
incident or the detected mode is not evident from Eq. 
(^) [it follows, however, from Eq. (371)]. 

The main conclusion from Eq. (39) is that one can- 
not neglect the matrices 012, 021. That they appear here 
demonstrates that the reflection and transmission prob- 
lem for > 1 are not directly related. It is tempting to 
interpret the additional fluctuations from these matrices 
as the origin of the quantity Z in Eq. (28). 

In the next subsection we discuss an intensity- weighted 
combination of all transmission delay times that docs not 
depend on 012 and 021. 



D. Weighted delay time and interpretation of long 
reflection delay times 

The matrix O, Eq. (^^, carries the correlations of 
the transmission channels and the eigenvectors of the 
Wigner-Smith matrix ('delay-time channels'). A suitable 
object which captures the essence of these correlations 
can be formed with help of the intensity-weighted delay 
times 



Im (dtnni/duj)t*^^ 
^ Ulm , I |2 



T ^1" * , 



2 

m I 1 



(40) 



where the last equality holds in the localized regime. The 
sum of all weighted delay times can be written as 



W = Im 



tr {dt/duj)^^ 



From Eq. 



tTtH 

we flnd the representation 



(41) 



+ 




2 3 
W/N(N+1)J 

FIG. 4. Distribution of intensity-weighted combination W 
of all transmission delay time for multi-channel wave guides 
with iV = 2 (full dots), N = 5 (open dots), and iV = 10 
(squares), from the numerical simulation. The full curve is 
distribution P(VFmax), Eq. (|^, of the upper bound Wmax- 
The dashed curve is the result for A'" = 10 if the matrices on 
and 022 in Eq. (^) would be random orthogonal matrices. 



W : 



tTT{oiiTO{i + 022r'022) 

2trr 

-[Oiifofi -I- 022t'o^2]i1' 



(42) 



where the first diagonal element is picked out because the 
transmission eigenvalue Ti is much larger than the other 
transmission eigenvalues. Hence W indeed carries infor- 
mation of the correlations between the dominant trans- 
mission channel and the delay-time channels, which can 
be quantified by the overlaps 



Oi 



On H, 



022 H- 



(43) 



Note that W does not involve the off-diagonal blocks 012 
and 021 of O which couple both ends of the wave guide, 
and that W is manifestly positive. 

The distribution P{W) is plotted in Fig. ^ for some 
values of A^ in units 'yN{N +1). These distributions are 
close to the distribution 



P(VF„,ax) 



[A(A-H)7]- 
VF3 

* ^ max 



■ exp 



Wn 



A(A + 1)7 



of the mean 



Winax = -(n 



(44) 



(45) 



of the two largest delay times fi, f[, which follows from 
Eq. (^. Fig. ^ also shows the distribution function if 
On and 022 would be random orthogonal matrices, which 
would result in much smaller values W ~ 7A^. 

The quantity Wmax is an upper bound of W. That the 
distributions of both quantities are very close requires a 
large overlap. 



7 



Ol 



1, 



(46) 



of the dominant transmission channel with the channel 
with the largest delay time, hence, that both channels 
are strongly correlated. The correlation appears to be 
most pronounced especially for large Wmax, because the 
tails of the two distributions coincide very well. 

Upon reflection, the strong correlations of the dom- 
inant transmission channel and the channel with the 
largest delay time can be seen as one reason why the 
single- mode delay times T„m are of order 'yN{N + 1), 
which corresponds to Y ~ Z ~ -/N{N +l)/2 in Eq. 



V. CONCLUDING REMARKS 

In this paper we have investigated the statistical prop- 
erties of the transmission delay time t in the presence 
of wave localization. Most of the analysis relied on the 
separation of the delay time into two independent con- 
tributions, T = X + x' 7 with X and x' given in Eq. (23). 
The properties of the delay time follow then from the 
distribution function of x and x' ■ This distribution does 
not depend on length in the localized regime. It is given 
as an exact analytic expression for = 1 in (|2l|) and in 
approximate form for A^ > 1 in Eq. ( p9|). 

We also have demonstrated in Sec. LVD that the dom- 



inant transmission channel is closely related with the 
channel associated to the largest Wigner-Smith delay 
time. Large reflection delay times can hence be inter- 
preted as exploration of regions deep inside the wave 
guide, which are only accessible via the dominant trans- 
mission channel. 

The separation rule (|2^) entails strict correlations 
among the delay times of a single realization, which are 
related by Eq. ( |24| ) . These relations become invalid when 
absorption dominates over localization (then diffusion 
theory becomes applicable again). It would be interest- 
ing to investigate whether the departure from Eq. ( |2^ ) 
qualifies as a practical tool that distinguishes these two 
distinct mechanisms of wave attenuation. 
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In this Appendix we briefly sketch how the evolution 
equation is be derived within one-dimensional scaling 
theory ||^,|2^,^,Q, adapted to the dynamical problem 
along the lines of Ref. |23y35| , |36|] . 

In this approach we study the evolution of x{L + 6L) = 
x{L) + 5x as the length of the wave guide is increased 
gradually, by adding a thin slice of length 5L. Within 
an ensemble of random disorder, the evolution of the dis- 
tribution function is then governed by a Fokker-Planck 
equation. 



(Al) 



In order to show that Eq. (Al) becomes Eq. (20), it 
remains to calculate the moments {5x) and {5x^). The 
scattering matrix elements 

ri = -ri* = iB, ti=t'^ = l + iA-{a + h) /2, (A2) 

of the slice are given by a Gaussian real number A with 
variance {A^) — a and the complex number B with 
(|B|2) = b. From (\r^\) = SL/21 we obtain the rela- 
tion to the mean free path / = dL/2b. The derivative 
dA/duj — 5L/c, as appropriate for the quasi-ballistic mo- 
tion through the small segment. 

Now we have to determine the elements u and du / duo 
for the composed system of length L + SL. From the 
composition rule (25) and the reduced form t = uvVt, 
Eq. (^), we obtain in the localized regime the prescription 

u{L + 5L) = u{l + iA + iRcBu^ - a/2 - 6/2), (A3) 

du , ^ du , , T , 

— (L + 6L) = —(1 + iA + iReBu^ - a/2-6/2) 

dio du 

, SL „ „ du\ 

u ( i h 2iReBu— , 

c duj J 



(A4) 



where we denoted for simplicity the initial value u{L) 
u. The increment of Sx is then given by 



c- SL ^ du 

Sx= h 2ReBu— , 

c du 



and the moments are 
SL 



{Sx) = -, {Sx')^2bx'^^-. 
c 7 c 



(A5) 



(A6) 



APPENDIX A: DISTRIBUTION FUNCTION OF x 

FOR N = 1 

The distribution P{x) ior N = 1, Eq. (|l|), follows from 
the requirement of stationarity on its evolution equation 
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